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Summary 

A method is proposed (without experimental verification) for 
extending the total strain version of Strainrange Partitioning 
(TS-SRP) to predict the lives of thermomechanical fatigue 
(TMF) cycles. The principal feature of TS-SRP is the 
determination of the time-temperature-waveshape dependent 
elastic strainrange versus life lines that are added subsequently 
to the classical inelastic strainrange versus life lines to form 
the total strainrange versus life relations. The procedure is 
based on a derived relation between “failure” and “flow” 
behavior. “Failure” behavior is represented by conventional 
SRP inelastic strainrange versus cyclic life relations, while 
“flow” behavior is captured in terms of the cyclic stress-strain 
response characteristics. Stress-strain response is calculated 
from simple equations developed from approximations to more 
complex cyclic constitutive models. For application to TMF 
life prediction, a new testing technique, bithermal cycling, is 
proposed as a means for generating the inelastic strainrange 
versus life relations. Flow relations for use in predicting TMF 
lives would normally be obtained from approximations to 
complex thermomechanical constitutive models. Bithermal 
flow testing is also proposed as an alternative to thermo- 
mechanical flow testing at low strainranges where the 
hysteresis loop is difficult to analyze. For the current paper, 
cyclic stress-strain response is calculated for the alloy 
2-l/4Cr-lMo steel in the post-weld, heat-treated condition for 
both thermomechanical and bithermal cycles and for a variety 
of strainranges and hold times. Bithermal failure behavior data 
are currently unavailable. 

Introduction 

Strainrange Partitioning (SRP) (ref. 1) was originally 
formulated on an inelastic strainrange versus cyclic life basis 
for isothermal conditions. The approach has worked well in 
the high-strain, low-life regime where the inelastic strains are 
large enough to be determined accurately by analytical and 
experimental methods. 

To extend the method into the low-strain, long-life regime 
where the inelastic strains are small and difficult to determine, 
it was necessary to consider the total rather than just the 
inelastic strainrange. Formulating SRP on a total strain basis 
(TS-SRP) (refs. 2 and 3) required the determination of both 


the inelastic and the elastic strainrange versus cyclic life 
relations for cycles involving creep strain. The elastic 
strainrange versus life relations for cycles involving creep are 
influenced significantly by temperature, hold time, waveshape 
(CC, CP, and PC), and how creep is introduced into the cycle 
(stress hold, strain hold, slow strain rate, etc.). To make the 
analysis tractable, we have assumed that the elastic life lines 
are displaced parallel to themselves and are parallel to the 
elastic life line for PP cycling. This means that for a given 
imposed cycle type at a specific temperature, the elastic line 
intercept (elastic strainrange at Nf - 1 cycle) is a function 
only of the time of the cycle. In the first version of TS-SRP 
(ref. 2), the elastic line intercept was determined by the use 
of an empirical equation with constants determined using data 
obtained from failure tests. Efforts to reduce testing require- 
ments led to the development of an updated version of TS-SRP 
(ref. 3). This development was based on a derived relation 
between failure behavior (inelastic strainrange versus cyclic 
life) and flow behavior (cyclic stress-strain response) and 
greatly reduced the time and cost of characterizing the TS-SRP 
behavior of an alloy. Failure behavior is characterized readily 
in the high-strain, low-life regime where testing times are 
reasonable, and the critical flow behavior is characterized in 
the desired low-strain, long-life regime by cycling a specimen 
just long enough for the stress-strain hysteresis loop to 
approach stability. With this approach, the elastic line intercept 
can be established with a minimum of long-time testing. 

Traditionally, thermomechanical fatigue 1 (TMF) resistance 
of materials has been estimated by conducting isothermal 
fatigue tests at the expected maximum temperature in the TMF 
cycle, with effects of creep being determined by imposing hold 
periods at peak tensile or compressive strain. Computer- 
controlled testing equipment and techniques have recently been 
developed to the point where thermomechanical fatigue tests 
can be conducted routinely, so isothermal testing is no longer 
as justifiable. However, true TMF cycles become difficult to 
control, analyze, and interpret at low strainranges, and 
bithermal fatigue testing has been proposed (ref. 4) to avoid 
this difficulty. 


^The term “thermomechanical fatigue” is used to indicate variable 
temperature fatigue in which the mechanical strain is imposed only by 
externally applied loads. Temperature gradients within the test volume are 
not permitted. 
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Existing methods for predicting TMF life of components 
in the low-strain regime are lacking in accuracy or generality. 
Thus, we have actively pursued development of a TMF life- 
prediction method that would overcome these deficiencies. As 
a starting point, we have selected the SRP method because 
of past successes. For example, the inelastic strain version of 
SRP has been successfully applied to a thermomechanical 
fatigue problem (ref. 5) at a strain level where the inelastic 
strainrange could be determined accurately and for which the 
life relations are temperature independent. Also, the TS-SRP 
version has been applied successfully to isothermal problems 
(refs. 2 and 3) where the inelastic strainrange cannot be 
determined accurately. These successes suggest that the total 
strain version has the potential to be extended to predict the 
lives of TMF cycles at low strainranges. 

This report examines the extension of TS-SRP into the TMF 
regime by the use of a hypothetical example for the pressure 
vessel and piping steel alloy 2-l/4Cr-lMo. The cyclic stress- 
strain-time, that is, flow response, has been determined for 
both thermomechanical and bithermal cycles using the 
Robinson flow model (ref. 6). Unfortunately, adequate TMF 
failure data for this alloy, especially in the low strain regime, 
are not available so we cannot verify the TS-SRP life predic- 
tions at this time. The characterization of an alloy by TS-SRP 
requires the determination of both failure and flow responses. 


Symbols 

A general constant in empirical flow equations 

A' general constant in empirical flow equations 

B intercept of elastic strainrange versus life relation 
C intercept of inelastic strainrange versus life relation 

C ' intercept of equivalent inelastic line for combined 

creep-fatigue cycles with parallel inelastic failure 
lines 

C" intercept of equivalent inelastic line for combined 
creep-fatigue cycles with nonparallel inelastic failure 
lines 

CC creep strain in tension, creep strain in compression 

CP creep strain in tension, plastic strain in compression 

/ frequency 

F strain fraction 

K cyclic strain-hardening coefficient 

N number of cycles to failure 

PC plastic strain in tension, creep strain in compression 

PP plastic strain in tension, plastic strain in compression 

R mean stress correction term for nonisothermal fatigue 

r correlation coefficient 

t hold time per cycle 


V fatigue mean stress correction term; isothermal and 
nonisothermal 

y general dependent variable 

A range of variable 

e strain 

e strainrate 

a stress 

r period of one cycle 

Subscripts: 

amp amplitude 

bi bithermal 

c compression 

cc creep strain in tension, creep strain in compression 
cp creep strain in tension, plastic strain in compression 

el elastic 

eff effective 

fm failure; mean stress condition 

f 0 failure; zero mean stress condition 

in inelastic 

ij pp, cc, pc, cp 

m mean 

max maximum value 

min minimum value 

pc plastic strain in tension, creep strain in compression 

pp plastic strain in tension, plastic strain in compression 

pre predicted 

t tension or total strain 

tm thermomechanical 

a stress 

Superscripts: 

b time-independent power of cyclic life for elastic 
strainrange versus life relation 

b' time-dependent power of cyclic life for elastic 
strainrange versus life relation 

c power of cyclic life for inelastic strainrange versus 

life relations 

c* power of cyclic life for inelastic strainrange versus 
life relations 

m general power of time in empirical flow equations 
n time-independent, cyclic strain-hardening exponent 

n f time-dependent, cyclic strain-hardening experiment 

a power on total strainrange in empirical flow equations 
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Cycle type (nonisothermal): 


HRIP high-rate inphase cycle producing 100-percent PP 
strainrange 

HROP high-rate out-of-phase cycle producing 100-percent 
PP strainrange 

THIP tensile hold inphase cycle with hold period at 
maximum tensile strain producing CP and PP 
strainranges 

CHOP compressive hold out-of-phase cycle with hold period 
at minimum compressive strain producing PC and 
PP strainranges 

Analysis 

In formulating TS-SRP for isothermal fatigue (refs. 2 and 3), 
we assumed that the inelastic and elastic strainrange versus 
life lines for creep-fatigue cycles are parallel to the 
corresponding lines for pure fatigue or PP cycles, as shown 
in figure 1. This is not an arbitrary assumption since our 
experience with several alloys suggests that this is reasonable 
behavior for isothermal conditions. Note that the frequency 
modified approach of Coffin (ref. 7) also uses parallel inelastic 
and elastic lines. We assume, in the absence of data for 
guidance, that this will also be the case for nonisothermal 
conditions. 

Based on the above assumptions, a relationship between 
failure behavior and flow behavior can be established. Failure 
behavior is expressed by the equations for elastic and inelastic 
strainrange versus cyclic life: 



0 


L0G(N f ) 



Figure 1 .—Relation between total strainrange and life for nonisothermal creep- 
fatigue cycles. Inelastic line intercept C' is determined from equation (3) 
and elastic line intercept B is determined from equation (7). 


explicitly consider the time dependency of the inelastic lines, 
although it could be added if needed. 

The four generic SRP relations are 

Ae in = Q/N'jY (4) 

The IDR is written as follows: 



~c> 

/ N 

II 

UJ 

<1 

(1) 

Ae in = C'(N fo ) c 

(2) 


where 

C' = [EFy(Q) l/c ] c (3) 

and 

ij = pp, cc, pc, or cp 

Equation (3) is derived from the interaction damage rule 
(IDR) of reference 8 and the four generic SRP inelastic 
strainrange— cyclic life relations for a theoretical zero mean 
stress condition. In the past the inelastic line intercepts for 
creep cycles (C pc , C cpy and C cc ) were taken to be independent 
of time. But recent developments (ref. 9) indicate that they 
may be time dependent at elevated temperatures, and 
procedures have been proposed for expressing the time 
dependencies analytically. The following development of the 
TMF life-prediction method based on TS-SRP will not 


where EFy = 1.0. Using equation (4) to solve for N t j and 
substituting into equation (5), we obtain equation (3). Flow 
behavior is expressed by an equation relating elastic and 
inelastic strainranges: 

Ae e/ = K i j(Ae in ) n (6) 

where n = b/c. 

Based on the assumption that the inelastic and elastic failure 
lines for creep-fatigue cycles are parallel to the corresponding 
failure lines for PP cycles, it follows that the strain-hardening 
exponent n in equation (6) is a constant as shown in figure 2. 

For isothermal conditions the strain-hardening coefficient 
K t j is a function of temperature, hold time, how creep is 
introduced into the cycle (stress hold, strain hold, etc.), and 
the strainrate-hardening characteristics of the alloy. For 
nonisothermal conditions it is also a function of the maximum 
and minimum temperatures and the phase relation between 
strain and temperature. 

The time-dependent behavior of the elastic strainrange-life 
relation for creep cycles is shown schematically in figure 3. 
Setting equation (1) equal to equation (6) and eliminating N fo 
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Figure 2.— Relation between inelastic and elastic strainranges for nonisothermal 
creep-fatigue cycles. Cyclic strain-hardening coefficient K )j is a function of 
hold time, maximum and minimum cycle temperatures, phase relation 
between strain and temperature, and strainrate-hardening characteristics of 
an alloy. 



Figure 3.— Relation between elastic strainrange and life for nonisothermal 
creep-fatigue cycles. Variation in elastic line intercept is a function of hold 
time, maximum and minimum cycle temperatures, phase relation between 
strain and temperature, and strainrate-hardening characteristics of an alloy. 


At r Ae^/ + Ae Jrt 


From equations (1) and (2) we obtain 

Ae, = B(N fo ) h + C'(N fo Y (9) 

A schematic plot of equation (9) is shown in figure 1 . Note 
that the solution of this equation gives the cyclic life Nf 0 for 
a theoretical zero mean stress condition. The final step in a 
life prediction is to adjust the computed life to account for any 
mean stress effects that may be present. 

A method for accounting for mean stress effects on life for 
isothermal conditions has been proposed in reference 10; 

(Nfrn) b = (N f0 ) b ~ V eff (10) 

where and Nj 0 are the lives with and without mean stress, 
respectively, and V eff is the effective mean stress correction 
term. For isothermal fatigue, V eff is determined by the 
following equation: 


Veff = Vo exp 



( 11 ) 


where V a — a m ! a amp . Note that this method of correcting for 
the mean stress effect on cyclic life was developed for a 
specific nickel-base alloy and may not apply to other alloys 
or even to other nickel-base alloys. The analyst must determine 
that the mean stress correction equation used is applicable to 
the alloy of interest. 

For TMF an alternate definition of V eff is in order since a 
mean stress can naturally develop because of the temperature 
dependency of the yield strength in tension and compression. 
Hence, V e ^ in equation (10) should be determined (ref. 11) 
by the following: 


using equation (2), we obtain the following equation relating 
flow and failure behavior: 

B = K i} {C’) n (7) 

In this equation the inelastic line intercepts C,y and the 
exponent c used to determine C' are considered to be failure 
terms. The strain fractions the strength coefficient K iJy 
and the strain-hardening exponent n are considered to be flow 
terms. Thus, the elastic line intercept B can be determined for 
a creep cycle from a combination of flow and failure data. 
Note that K tj and F tj will, in general, depend on waveform. 

We are now in a position to establish a total strainrange 
versus life relation and thus predict life on a total strainrange 
basis. Note that the SRP inelastic strainrange versus life 
relations and the flow relations are determined for a specified 
minimum and maximum temperature and phase relationship 
between strain and temperature. The total strainrange is 


V eff = — (12) 

1 -* 

R y 

where R a is equal to a min /a max and R v is the absolute value of 
the ratio of the compressive yield strength to the tensile yield 
strength at their respective maximum and minimum 
temperatures and strainrates in the TMF cycle. To date there 
is no direct experimental verification of this method for 
accounting for mean stress effects for nonisothermal fatigue. 

To predict life on a total strainrange basis, it is necessary 
first to determine the PP inelastic and elastic lines and the 
desired SRP inelastic strainrange versus life relations 
experimentally. Note that these relations are to be established 
for a theoretical zero mean stress condition. Empirical 
estimation methods developed for isothermal fatigue, such as 
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the ductility-normalized SRP relations (ref. 12) are not 
recommended at this time as they have not been verified for 
application to TMF. Ideally, failure and flow behavior would 
be determined from TMF tests duplicating the cycles for which 
lives are to be predicted. However, this approach is impractical 
as it lacks generality of use. For example, if the cycle were 
to change, the entire data base would have to be regenerated 
at a doubling of cost and lead time. Further changes would 
in turn require further repetition of experiments. A more basic 
approach is thus required. While an isothermal approach would 
offer advantages in terms of costs because of the vast 
background of isothermal data bases, we do not, at this time, 
recommend doing so. This recommendation comes as a result 
of an indepth survey (ref. 11) comparing the TMF and 
isothermal fatigue resistances of many alloys. Only under 
special circumstances of temperature invariant deformation and 
cracking mechanisms could isothermal fatigue resistance be 
used to accurately predict TMF results. 

At this stage of development of the TS-SRP approach, we 
recommend resorting to bithermal fatigue tests (ref. 4) to 
generate the experimental inelastic strainrange-life relations 
required by TS-SRP. Bithermal fatigue offers the simplicity 
of isothermal testing, yet it captures the first-order effects of 
inphase and out-of-phase TMF cycling. The bithermal tests 
should cover a sufficient temperature range to encompass the 
deformation and cracking mechanisms pertinent to the TMF 
cycles of interest. Determination of the critical temperature 
ranges for testing will require a rudimentary understanding 
of the metallurgical factors governing the deformation and 
cracking mechanisms. Thermomechanical flow tests would 
normally be conducted to properly characterize the stress-strain 


response, but bithermal flow tests could be conducted at the 
lower strainranges where the inelastic strains are small and 
it is difficult to analyze the thermomechanical hysteresis loops. 
The stress- strain response of the two types of cycles should 
be similar in this strain regime, and the bithermal cycle would 
be a good approximation to the thermomechanical cycle. The 
thermomechanical and bithermal cycles used in this report are 
shown in figure 4. 

Techniques are described in reference 4 for determining PP 
life relations for inphase and out-of-phase bithermal cycles, 
CP inphase, and PC out-of-phase inelastic strainrange-life 
relations. Proposals for the determination of a CC bithermal 
life relation have not been discussed because of the exclusion 
of such a strainrange component in TMF cycles at small 
inelastic strainranges. 

The strain-hardening coefficient K tJ and the strain fraction 
Fjj can be determined using an appropriate constitutive flow 
model for which the material constants are available. As an 
alternative, they could be determined by conducting flow tests 
for the creep-fatigue cycles of interest. Using these data, one 
can determine the necessary empirical correlations describing 
flow behavior. This latter approach is the most viable option 
at this time because reliable constitutive flow models in the 
low strain regime and the required material constants are not 
yet generally available. The procedures for determining the 
flow correlations are described in the following section. 

The strainrange analysis presented above is based on the 
assumption that the inelastic and elastic versus life lines for 
creep-fatigue cycles are parallel to the corresponding lines for 
PP cycles. This may not always be a satisfactory assumption. 
(The case of nonparallel lines is discussed in the appendix.) 




(a) PP cycle, inphase. (b) PP cycle, out-of-phase, 

(c) PC cycle, out-of-phase. (d) CP cycle, inphase. 


Figure 4. —Bithermal and TMF waveshapes used in this report. 
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Analysis Using Robinson Constitutive 
Model 

The choice of a constitutive model for use with TS-SRP 
is somewhat arbitrary. We have selected the Robinson model 
(refs. 6 and 13) because it has been validated for TMF 
application (using the alloy 2-l/4Cr-lMo steel in the post- weld, 
heat-treated condition). Using Robinson’s model, we have 
obtained the simple power law correlation shown below. This 
same power law form was also used successfully to correlate 
isothermal flow data (ref. 3): 

y = A(t) m (13) 

where y is the dependent variable representing several different 
flow variables, as will be discussed shortly, and t is the hold 
time per cycle. 

Generally, the intercept A (value of y at t = 1) is a function 
of total strainrange. The results obtained from the Robinson 
model for thermomechanical cycles and earlier results for 
isothermal cycles using the Walker model (refs. 3 and 13) show 
that the family of lines shown in figure 5 can be taken as 
parallel. Thus, the exponent on time m is assumed to be 
independent of total strainrange. By a process of trial and error, 
we determined that the intercept A can be correlated with total 
strainrange by another power law as shown in figure 6; 



Figure 5.— Power law relation used to correlate flow data. Lines are parallel, 
and intercept A at t = 1 sec is a function of total strainrange. 



A=A'( Ae t ) a 

thus 

y = A'(Ae) a (t) m (15) 

The dependent variable y is now a function of two inde- 
pendent variables, Ae t and t. If both sides of equation (15) are 
divided by (Ae,) a , the family of lines shown schematically in 
figure 5 collapse to the single line of figure 7. The values of 
A', a, and m vary with the dependent variable y and the 
mechanical properties of the alloy. Note that TS-SRP is not 
dependent on the form of the equation used to correlate the 
flow data and that equation (15) could be of many different 
forms. The only requirement is that it represent the flow data 
in a tractable form with sufficient accuracy. 

Five flow correlations based on equation (15) are used herein 
to determine the required flow variables: K ij} Fy, A <r, o t , and 
o c . The first two are used to determine the coefficients B and 
C' in equation (9). The remaining three are used to determine 
the term R a in the mean stress correction using equation (12). 

Note that, in principle, each of these correlations could be 
obtained directly from a suitable constitutive model. Although 
the exact form of the relations would no doubt differ somewhat 
from the one used here, the trends would be quite similar. 
The empirical correlations are used only because of their 
extreme simplicity and comparatively good accuracy. 

For a specific alloy, these correlations depend on the 
maximum and minimum temperature of the cycle, the 
waveshape, how creep is introduced into the cycle (stress hold, 
strain hold, etc.), the straining rate during loading and 
unloading, and the phase relation between strain and 
temperature. Only inphase and out-of-phase continuous cycles 
and strain-hold cycles with zero mean strain are considered 
at the moment. These are shown in figure 4. When computing 
the stresses and strains for the postulated bithermal hysteresis 
loops, the temperature is changed at zero stress on the tension- 
going and compression-going sides of the loops. And when 
calculating the TMF loops, the temperature ramp rate is 
determined by the maximum and minimum temperatures in 
the cycle and the mechanical straining rate. The relation 
between strainrate and temperature rate for a sawtooth 



Figure 7.— Power law normalized on total strainrange raised to suitable power, 
collapsing family of lines shown in figure 5. 
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waveform is derived in the following manner for inphase and 
out-of-phase cycles. For a constant strainrate, 

e = 2/ ( Ae,) (16) 

The period of one cycle is the reciprocal of the frequency / 
(r = 1 If), and the temperature range is traversed in one-half 
the cycle period. 

Let t' equal r/2, thus the time required to go from one strain 
limit to the other is 


r' 



(17) 


Let AT equal the temperature range of the cycle of interest 
(inphase or out-of-phase). The temperature ramp rate T is given 
by the temperature range divided by the time required to 
traverse the cycle. Note that t is positive if the temperature 
is increasing and negative if the temperature is decreasing; 


^ A7IAe,) 

e 


(18) 


The inphase and out-of-phase thermomechanical and 
bithermal hysteresis loop results obtained from the Robinson 
model are given in table I and II, respectively, for the following 
conditions: 

Total strainrange values 

Continuous cycling: 0.002 to 0.010 in increments of 0.001 
Strain-hold cycling: 0.002, 0.003 , 0.004, 0.006, 0.008, 
0.010 

Hold time, sec: 60, 300, 600, 1800, 3600 
Minimum temperature, °C: 250 
Maximum temperature, °C: 600 
Strainrate, in. /in. /min: 0.04 

A review of the results of these computations for the cycles 
considered herein reveals the following differences in the 
stress-strain response of the thermomechanical and bithermal 
cycles for a given total strainrange and hold time. 


Inphase cycles 

^ \ a t,bi\ 

I ^c,tm 1 ^ \®c,bi \ 

|Aa, m | < \Ao bi \ 

I &€el ,tm I < |A€*/,wl 


Out-of-phase cycles 
^ Kbi\ 

\®c,bi I 

I Acr fm | < |Aa w | 

I I ^ |A^ el,bi I 


The above results are as expected. For example, in an 
inphase bithermal cycle, the stress is acting at the maximum 
temperature for the entire duration of the tensile half, but in 
a thermomechanical cycle the stress is acting at the maximum 
temperature only briefly at the maximum strain limit. Thus 


o t will be greater in a thermomechanical cycle than in a 
bithermal cycle because of the greater thermal recovery in the 
bithermal cycle. A similarly based argument also applies to o c . 

We are now able to determine the flow correlations listed 
previously using equation (15), but first the strain-hardening 
exponent n in equation (6) must be determined using time- 
independent PP flow data. The strainrate of 0.04/min may not 
be fast enough to obtain time-independent (PP) deformation 
during loading and unloading but was the fastest rate used when 
the material constants for the model were determined. We have 
assumed that this rate produces no time-dependent inelastic 
strains. Note also that it is highly unlikely that the high- 
temperature ramp rate T implied by equation (18) could be 
achieved during the cooling leg of a thermomechanical cycle. 
Nevertheless, the following results were obtained: 

A € e i = 0.0045(Ae, n )° 107 (19) 

for thermomechanical cycles and 

Ae el = 0.0044(A€,„)° 105 (20) 

for bithermal cycles. The relations for inphase and out-of-phase 
PP cycles are identical. 

The differences between the thermomechanical and bithermal 
results are virtually nill under the current circumstances. These 
correlations are shown in figure 8 wherein the symbols 
represent values calculated using the Robinson model and the 
line represents the empirical correlation. Correlation coeffi- 
cients of 0.990 and 0.995 for the thermomechanical and 
bithermal cycles, respectively, indicate an exceptionally good 
representation of the Robinson model. In the following, we 
present only results for TMF cycling. 

Correlation Between Cyclic Strain-Hardening Coefficient 
and Hold Time 

Equation (15) expressed for K t j is as follows: 


K^A'iAe^itr (21) 

The strain-hardening coefficient for each loading condition 
considered using the Robinson model (hold time, total 
strainrange, phase relation, etc.) is calculated using equation (6) 
and the proper value of the strain-hardening exponent from 
equations (19) and (20). Previously (ref. 3), K t j was taken to 
be independent of total strainrange, but additional analysis 
indicates that it is a weak function of total strainrange. 
A multiple regression analysis of the appropriate values gives 
the following correlations for out-of-phase and inphase TMF 
cycling with n = 0.107: 

_*£L_ = 4. 68 9 x 10 ~\t) -° 0167 (22) 

(Ae r ) 
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(a) Thermomechanical cycles; inphase and out-of-phase. 

(b) Bithermal cycles; inphase and out-of-phase. 

Figure 8.— Relation between elastic strainrange and inelastic strainrange for 
time-independent (PP) deformation for nonisothermal cycles. Alloy, 
2-l/4Cr-lMo steel in post-weld, heat-treated condition; minimum temper- 
ature, 250 °C, maximum temperature, 600 °C. 


(a) PC cycles, out-of-phase. 

(b) CP cycles, inphase. 

Figure 9.— Relation between strain-hardening coefficient and hold time 

for thermomechanical strain-hold cycles. Alloy, 2-l/4Cr-lMo steel in post- 
weld, heat-treated condition; minimum temperaure, 250 °C, maximum 
temperature, 600 °C; cyclic strain-hardening exponent, n , 0.107. 


for out-of-phase cycling and However, good correlations were obtained by dividing the data 

into two regimes as indicated below. 

K For 0.002 < Ae, < 0.004, 

= 5.052 x 10 - 3 W-° 0158 (23) 

(AC() 

f-2 073 = 1-416X 10 -6 (0° 0506 (25) 

for inphase cycling. These equations are represented in figure 9, (Ae,) 

and their excellent ability to correlate the calculated results 

is shown in figure 10. f or out-of-phase cycles and 

Correlation Between Strain Fraction and Hold Time 

F 

Equation (15) expressed for Fy is as follows: (Ae )~ 2110 = 10 6 (0 00448 (26) 

Fij “ A' (Ae t ) a (t) m Fij < 1.0 (24) 

for inphase cycles. These equations are expressed in figure 11, 
Analysis of the strain fraction— hold time data showed that and their ability to correlate the calculated results is represented 

equation (24) is not applicable over the entire range of the data. by correlation coefficients of 0.993 and 0.994, respectively. 
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(a) PC cycles, out-of-phase. 

(b) CP cycles, inphase. 

Figure 10.— Ability of equation (21) to correlate cyclic strain-hardening 
coefficient with hold time for thermomechanical cycles. Alloy, 
2-l/4Cr-lMo steel in post-weld, heat-treated condition; minimum temper- 
ature 250 °C; maximum temperature, 600 °C. 

For 0.004 < Ae t < 0.010, 


= 6.062 X 10 - 5 (t) 00744 (27) 

(Ae,) - 367 

for out-of-phase cycles and 

— = 6 ' 166X 10- 5 (0°° 733 (28) 


for inphase cycles. These equations are expressed in figure 12, 
and the correlation coefficients are 0.999 and 0.997, 
respectively. The constants for the correlations for K tj and F )j 
for TMF cycling (eqs. (21) to (28)) are summarized in table HI. 

Correlations Between Stress and Hold Time 

Our experience suggests that better correlations for a t and 
a c are usually obtained when o t is used for cycles where creep 



(a) PC cycles, out-of-phase. 

(b) CP cycles, inphase. 

Figure 11.— Relation between strain fraction and hold time for thermo- 
mechanical strain-hold cycles. Alloy, 2-l/4Cr-lMo steel in post- weld, heat- 
treated condition; minimum temperature, 250 °C, maximum temperature, 
600 °C; 0.002 < Ae, < 0.004. 

occurs in the tensile half of the hysteresis loop (CP cycle) and 
when o c is used where creep occurs on the compressive side 
of the loop (PC cycle). The results obtained here show that 
stress is a very weak function of hold time and could be omitted 
with little loss of accuracy. But this may not be true generally, 
and we have chosen to include it for illustrative purposes. The 
resulting stress correlations are summarized in table IV. 

Life Prediction of TMF 

In this section we outline the steps required to predict the 
life of a TMF cycle. For purposes of illustration, an inphase 
tensile strain-hold cycle (THIP) is selected. The mechanical 
strainrate and temperature limits are given in the previous 
section. 

TMF cycles invariably contain both time-dependent and 
time-independent components of inelastic strain. Thus a THIP 
cycle will contain both PP and CP strain components, and the 
appropriate generic SRP inelastic strainrange versus life 
relations, equation (4), are required. Both relations must be 
for inphase cycling. Unfortunately, the numerical values of 
the material constants Qj and c are not available at present. 


9 



HOLD TIME, SEC 

(a) PC cycles, out-of-phase. 

(b) CP cycles, inphase. 

Figure 12.— Relation between strain fraction F » and hold time for thermo- 
mechanical cycles. Alloy, 2-l/4Cr-lMo steel in post-weld, heat-treated 
condition; minimum temperature, 250 °C; maximum temperature, 600 °C; 
0.004 < Ac, < 0.010. 

As discussed earlier, bithermal testing is recommended for 
determination of these material constants. For now, we 
proceed as if Cy and c values are known; 


Ae in = C pp (N pp Y 

(29) 

A e in = C cp (N cp ) c 

(30) 


The intercept C' of the equivalent inelastic line in equation (9) 
can now be determined. From equation (3), 

C ' = | F pp (C pp ) llc + F cp (C cp ) ,/c ] c (31) 

Since F pp - LO - F cp , equation (3) can be rewritten as 
follows: 



The strain fraction F cp is determined using the appropriate 
correlation given in the previous section. 

If F cp is very small (-0), C' = C pp ; and if F cp approaches 
unity, C' - C cp . The equivalent elastic line intercept B in 
equation (9) can now be determined using equation (7). The 
value of K cp is determined using the correlation given in the 
previous section, and the value of C' is determined using 
equation (32); 


B = K cp (C') n (33) 

The ingredients required to make a life prediction are now 
available. With a knowledge of Ae, and the proper constants 
for equation (9), we can now solve for N fo . The value of Nj 0 
can be determined by trial-and-error or by direct use of the 
inversion method given in reference 14. Note that Np is the 
cyclic life for a theoretical zero mean stress condition (ref. 10). 

The final step in a life prediction is to account for the effects 
of mean stress on cyclic life. Rewriting equation (10) yields 

Nfin = [(N fo ) b - V eff ] vb (34) 

The value of F^can be determined using equation (12) or 
some other method. If equation (12) is used, the values of 
a m in and are obtained using the appropriate stress versus 

hold time correlations given by the constants in table IV. 


Concluding Remarks 

The total strainrange version of Strainrange Partitioning 
(TS-SRP) was developed originally for isothermal fatigue. 
This development makes it easier to characterize an alloy and 
predict cyclic life in the low-strain regime without having to 
conduct failure tests in this regime. This development is based 
on a derived relation between failure behavior and the cyclic 
stress-strain or flow response of an alloy. Failure testing is 
done only in the high-strain regime where test times and costs 
are more reasonable, and flow testing is done in both the high- 
and low-strain regime. The flow tests need be run only until 
the stress-strain hysteresis loop approaches cyclic stability. 
Failure tests should also be done to determine the effects of 
mean stress on cyclic life. If mean stress effects are not 
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accounted for, inaccurate life predictions will result from using 
this or any other life-prediction method in the low-strain, long- 
life regime. 

A method for extending TS-SRP to characterize an alloy 
and predict the lives of thermomechanical cycles in the low- 
strain regime is presented. This method is based on noniso- 
thermal rather than isothermal data. Bithermal fatigue testing 
is recommended at this time to generate the nonisothermal data 
required to determine the inelastic strainrange-life relations 
required by TS-SRP. Bithermal fatigue testing offers the 
simplicity of isothermal testing but captures the first-order 
effects of inphase and out-of-phase TMF cycling. Thermo- 
mechanical flow tests would normally be conducted to properly 
characterize the stress-strain response, but bithermal flow tests 
could be conducted at the lower strainranges where the inelastic 


strains are small, and it is difficult to analyze the thermo- 
mechanical hysteresis loops. 

The data from flow testing are used to determine flow 
response of an alloy. This response is approximated using a 
simple power law with two independent variables with 
constants determined by a multiple regression analysis. These 
same power law relations were used in the isothermal version 
of TS-SRP. It is not a requirement that this particular power 
law be used, since any relation could be used provided it 
accurately represents the data. Currently, there are no TMF 
data to validate the proposed life-prediction method. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 4, 1987 
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Appendix— Analysis for Nonparallel Failure Lines 


Our experience strongly suggests that the assumption that 
inelastic and elastic strainrange versus life lines for creep- 
fatigue cycles are parallel to the corresponding lines for pure 
fatigue or PP cycles is reasonable for isothermal creep-fatigue 
analysis. On this basis, we have extended this assumption to 
nonisothermal creep-fatigue. However, this assumption may 
not be applicable for all alloys, and a method would be required 
for dealing with such a problem. 

If the assumption of parallel lines is not reasonable, then 
the equations relating failure and flow behavior given here are 
not applicable. For purposes of discussion, let us assume that 
the inelastic failure lines for creep-fatigue cycles (PC, CP, 
CC) cycles are not parallel to the PP line or to each other and 

that the elastic failure lines for these cycles both translate and 

rotate as a function of hold time. 

Failure behavior is now expressed by the following 
equations: 

Ae e ; = B(N fo ) y (35) 

where 

b' =b'(t ) 

The four generic SRP life relations are 

Ae* = C pp (N pp Y (36) 

for PP cycles and 


and solving for Ae in we obtain 

A e in = C"(N fo ) c (40) 

where 

c* = [F pp (C pp ) Uc + Fij(Cij) 1/c '(Ae,„) 1/c ~ 1/f * (41) 

Flow behavior is expressed by the following equation: 

Ae^^AeJ"' (43) 

Setting equation (35) equal to equation (43) and eliminating 
Nf 0 using equation (40) we obtain 

B = K ij (C") n ' (44) 

where n' = b' /c. 

Since the total strain version of Strainrange Partitioning 
(TS-SRP) is intended for the low-strain regime where the 
elastic strainrange is the dominate term in equation (8), it is 
reasonable to retain the assumption of parallel inelastic lines 
and to allow the elastic line to both translate and rotate with 
hold time. Assuming the inelastic lines to be parallel simplifies 
the analysis considerably, but does not introduce large errors. 
Thus equation (9) becomes 

Af, = B(N fo ) b + C'(N fo ) c (45) 


Ae,„ = QjM/' (37) 

for PC, CP, and CC cycles, where c* is not equal to c and 
c* may be different for each type of cycle. 

The interaction damage rule (IDR) is written as follows for 
a creep-fatigue cycle: 


Fpp j Fjj _ 1 

Npp ^ ij N fo 


(38) 


Using equations (36) and (37), we obtain the following: 


1 pp 




_ 1 _ 
A Tf 0 


(39) 


where b ' = b ' (f). Since b ' is time dependent, it then foil vs 
that the strain-hardening exponent n ' in equation (43) is also 
time dependent. The Ae d - Ae in line both translates and 
rotates with time. Thus, equation (15) cannot be used to 
determine K {j since this equation accounts only for 
translation. 

The time dependency of the fatigue exponent b ' in equation 
(45) can be determined by conducting flow tests to characterize 
the time dependency of the strain-hardening exponent n ' in 
equation (43). At this time we are unable to suggest a suitable 
form for the equations for K t j and n ' , but we feel that 
equation (15) can still be used to correlate the strain fractions 
Fij and the stress terms Aa, a„ and o c . 

When the life relations are nonparallel, failure testing will 
have to include tests at lower strainranges to characterize this 
more complex material behavior. 
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TABLE I. — RESULTS OF FLOW CALCULATIONS FOR THERMOMECHANICAL CYCLES USING THE ROBINSON MODEL 
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TABLE III. — CONSTANTS FOR K i} AND F tj CORRELATIONS FOR THERMO- 
MECHANICAL STRAIN-HOLD CYCLING FOR y = A '(Ae,) a (/) m 


[Material, 2-l/4Cr-lMo steel; post-weld, heat-treated condition.] 


Cycle type 

Total 

strain- 

range, 

Ae, 

Flow 

variable, 

y 

Constant, 

A' 

Power on 
total 
strain- 
range, 
a 

Power 
of time, 
m 

Correlation 

coefficient, 

r 

PC 


K P c 

4.689X 10 -3 

0.020 

-0.0167 

0.857 

(out-of-phase) 







0.002 
to 0.004 

F pc 

1.416X10" 6 

-2.073 

.0506 

.993 


0.004 
to 0.010 

F pc 

6.062 X10 -5 

-1.367 

.0744 

.999 

CP 


Kc P 

5.052 xlO" 3 

0.037 

-0.0158 

0.929 

(inphase) 







0.002 
to 0.004 

Fc P 

1.204X10 -6 

-2.110 

.0448 

.994 


: 0.004 
to 0.010 

F cp 

6.166X10- 5 

-1.364 

.0733 

.997 


TABLE IV. — CONSTANTS 3 FOR STRESS CORRELATIONS FOR 
THERMOMECHANICAL STRAIN-HOLD CYCLING FOR 
y = A'(Ae t nt) m 

[Material, 2-l/4Cr-lMo steel; post- weld, heat-treated condition.] 


Cycle 

type 

Flow 

variable, 

y 

Constant, 

A* 

Power on 
total 
strain- 
range, 

OL 

Power 
of time, 
m 

Correlation 

coefficient, 

r 

PC 

A o 

1394.7 

0.230 

0.0008 

0.990 


°c 

722.3 

.183 

.0019 

.989 

CP 

Act 

1379.8 

0.227 

0.0008 

0.987 


°t 

781.4 

.329 

.0007 

.986 


a Stresses in units of MPa (1 ksi = 6,895 MPa). 
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